Compactons in PT-symmetric generalized Korteweg-de Vries Equations 
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In an earlier paper Cooper, Shepard, and Sodano introduced a generalized KdV equation that 
can exhibit the kinds of compacton solitary waves that were first seen in equations studied by 
Rosenau and Hyman. This paper considers the PT-symmetric extensions of the equations examined 
by Cooper, Shepard, and Sodano. From the scaling properties of the "PT-symmetric equations a 
general theorem relating the energy, momentum, and velocity of any solitary-wave solution of the 
generalized KdV equation is derived, and it is shown that the velocity of the solitons is determined 
by their amplitude, width, and momentum. 



I. INTRODUCTION 

In a previous investigation Cooper, Shepard, and 
Sodanoi introduced and Khare and Cooper— studied fur- 
ther the first-order Lagrangian 



L{l,p) 



dx 



(1) 



This Lagrangian is a modification of the original com- 
pacton equations of Rosenau and Hyman^,. (Unless oth- 
erwise specified, the range of x integration is over the 
entire real axis — c» < a: < oo.) This Lagrangian gives 
rise to a general class of KdV equations of the form 



Ut 



Apu- 



(2) 



where the solution u{x,t) to the generalized KdV equa- 
tion is defined by u{x,t) = (px{x,t). For < p < 2 
and I ~ p + 2 these models admit compacton solutions 
whose width is independent of the amplitude. For p > 2 
the derivatives of the solution are not finite at the bound- 
aries of the compacton where u —^ 0. Cooper, Khare, and 
Saxenai^ analyzed the stability of the general compacton 
solutions of this equation and showed that solutions are 
stable provided that 



2 < ; <p-|-6. 



(3) 



This extension of the KdV equation is not a Hamiltonian 
dynamical system at arbitrary e. A more recent study 
by Fring^ was based on a Hamiltonian formulation. The 
Hamiltonian studied by Fring is related to a special case 
of the system of generalized KdV equations examined 
here. 

To find extensions of the generalized KdV equation 
that are invariant under the joint operation of space re- 
flection (parity) V and time reversal T, we make the fol- 
lowing definitions: spatial reflection V consists of making 
the replacement x —f —x. Also, because u is a velocity, 
under V we replace u by ~u. The effect of the time re- 
versal operation T is to change the signs of i, t, and u: 
i —i, t —t, and u ~u. Therefore, the combina- 
tion iux is VT even, so a T'T-symmetric generalization 
of the Lagrangian ([T|) is 



L-pT ^ dx 



1(1^1) 



(5) 



For this Lagrangian we must find the correct PT- 
symmetric contour that lies on the real axis when m — 2. 
For VT to be a good symmetry, branch cuts must be 
taken along the positive imaginary axis in the complex- a; 
plane. The Hamiltonian resulting from the above La- 
grangian is 



H 



dx 



1(1-1) 



avF{iUxy 



(6) 



There has been some recent interest in complex VT- 
symmetric extensions of the ordinary KdV equation. 
Such extensions exist in the complex plane but also lead 
to new PDEs that are entirely real. The first extension 
of the KdV equation by Bender et al^ was 



Ut - iu{iUxY 



0, 



which reduces to the usual KdV equation when e : 
This equation was analyzed by Bender et al^ for e 



(4) 

= 1. 
= 3. 



When m is an even integer, a convenient choice for a 
that allows for solitary-wave solutions and that gives a 
real equation for the generalized KdV system is 



a{m — l)i'' 



1. 



(7) 



For simplicity, we choose a as in ([7]) for most of this pa- 
per. The VT generalization of ([2]) has the same canonical 
structure as the KdV equation. From Lagrange's equa- 
tions or from Hamilton's equations for the generalized 



2 



KdV equations, we obtain the equations of motion for 
t): 



du d SH 



{u,H}, 



dt dx Su 
where the Poisson bracket structure iai 

{u{x),u{y)} = dj{x ~ y). 

The resuhing equation becomes 



(8) 



(9) 



0^ut+ u.^u^-"^ + uP-2ym-3 „ 2)r 
+2mpuuxxul + mu'^UxxxUx + {p - ■ (10) 

This system of equations has three obvious conserva- 
tion laws: conservation of mass M, momentum P, and 
energy E, where the energy is the value of the Hamilto- 
nian ^ and 



M 



dxu{x,t), P 



dx^u^{x,t). (11) 



The case m = 2 leads to the well known compacton so- 
lutions. Fring^ studied a Hamiltonian similar to the sub- 
class of this PT-symmetric class of Hamiltonians corre- 
sponding to / — 3 and p = 0, but with slightly different 
coefficients for the two terms in the Hamiltonian. 

This paper is organized as follows: Section |TT] consid- 
ers the scaling properties of the nonlinear wave equation 
(|10p and discusses the energy and momentum of solitary 
waves, and Sec. IIIII continues the discussion of these con- 
served quantities. Traveling-wave solutions are discussed 
in Sec. IIVI and the conserved quantities for these solu- 
tions are examined in Sec. |Vl Some special cases are de- 
scribed in Sec. IVII The question of stability is addressed 
in Sees. EIIl and IVml 



II. SCALING PROPERTIES 

Let us examine the scaling properties of pH)) . We re- 
quire that solutions transform into solutions under the 
scaling 

x^ Xx, X'H, u X^u, (12) 

and we find that 

/3 - r/ = (? - l)/3 - 1 = /3(p + m - 1) - m - 1. (13) 
Solving for (3 we obtain 



p + m — I 



We also find that 



1 - 7? = p{l - 2). 



(14) 



(15) 



Suppose that we have a traveling solitary wave of the 
ioim f{x—ct). Then c scales as x/f or as A^^**. Therefore, 
c scales as 

In terms of the velocity, x scales as A and 

p + m — I 



(X c , Zi 



m{l - 2) 



(17) 



From the equation for ii we see that the width of the 
solitary wave does not depend on the velocity when 



I = p + m. 



(18) 



This is a generalization of the result for to = 2. The 
conserved momentum scales like u'^x oc SO 



(X C"^ 7,2 = 



3r7i — I + p 



(19) 



m{l-2) 

Finally, the conserved energy scales as u'^x, so 

E cx A'''+^ cx c('"'+P+'"-'))/(p+'"-0, (20) 
Eliminating c in favor of P in this formula, we get 

EoiP-\ (21) 

where 

Im + p + m — I 



p + 3r77 — I 



(22) 



This reduces to our previous results for the case to = 2. 

We show below how to make these scaling laws more 
precise and how to determine the constants of propor- 
tionality that relate the conserved quantities. 



III. ANOTHER WAY TO RELATE ENERGY 
AND MOMENTUM OF SOLITARY WAVES 

In Ref. [4] a general theorem was derived that relates 
the energy, momentum, and velocity of solitary waves of 
the generic form u{x,t) = AZ[(3{x — q{t))] for to = 2. 
Here, we generalize this result to arbitrary to and we use 
the results obtained by studying the scaling properties of 
the alternative action $ to justify the previous derivation. 

We start from the action 



dtL, 



(23) 



where L is given in ([5]). We now assume that the exact 
solitary-wave solution has the generic form 



cf>x = u = AZ[P{x ~ ct)]. 



(24) 



Using this form, it is easy to calculate the value of the 
Hamiltonian ^ for the solitary wave 



H 



-Ci{l) 



I3l{l - 1) 



+ Af+™/3™-^C2(p,TO), (25) 
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dzZ\z), C2(p,m)= dz[Z'{z)]"'ZP 



where 

(26) 

Since H and momentum P are conserved, we can rewrite 
the parameter A in terms of P: 



IV. TRAVELING- WAVE SOLUTIONS 



1 r a2 

P=- dxu^ = — C5, 

2 J 2/3 ^' 



y" dz z^{z) 



(27) 



Replacing A by P, we rewrite the Hamiltonian H as 



(28) 



where 



C3(0 

C4(p, m) 



aC,{p,m) (2/C5)(^+'")/^ 



(29) 



At this point, we note that the exact solutions have the 
property that they are the functions of the parameter /3 
that minimize the Hamiltonian with respect to P when 
the momentum P is fixed. Using dH/d(3 ~ 0, we obtain 



This leads to 



C4{p, m){p + 3m — 2) 



where r is given by 



H^f{l,p,m)p-'-, 
and 

p + m — I 



(30) 



(31) 



p + 3rn — 2 
C4(p, m)(p + 3m — 2) 



(32) 



C3{l){l-2) 

Hamilton's equation q = dH/dP yields the relationship 

-q = c = rH/P. (33) 

From this analysis it is again easy to show that the mo- 
mentum P, amplitude and width parameter /? func- 
tionally depend on the velocity c (note that c = —q)'- 



p + 3m-l 1 l-p-m 



(34) 



Here, the proportionality constants depend on Ci{l,p, m) 
(i = 1, 2, 3, 4, 5) defined above, and once an exact solution 
is obtained, these constants can be calculated easily. 

We make three observations: (i) When I ~ p + m, 
the width parameter f3 is independent of the velocity c 
and momentum P and hence of the amplitude A of the 
solitary wave, (ii) The c dependence of the amplitude A 
depends solely on the parameter I, and it is independent 
of the parameters p and m. (iii) The stability problem 
when I ^ p + 2 was studied by Dey and Khare^ for the 
case m = 2 using the results of Karpman^^. In Sec. VII 
we extend their results to arbitrary even integer m. 



We begin with the wave equation pO|) which can be 
reexpressed as 



m — 1 '■ 



(uf<-i) = 0, 



m — 1 
and assume that 

u{x,t) = f{x - ct) = f{y). 

Then, 

cf = f-'f + 



(35) 



(36) 



{p [F'HfT]'- m [nrr^-']") , (3?) 



and integrating once we obtain 



cf = i—+mr{fr-'f" +p{f'rF~'+i<u m) 

where /' = df /dy. 

For compact solutions Ki = 0. Setting the integration 
constant Ki to zero, multiplying this equation by /, and 
integrating over y, we obtain 



1 p + m 
ch = -// 7 J, 

/ — 1 m — 1 



(39) 



where 

poo 

in= dyriy), j^,p= I dy[fTr{y)- m 



For noncompact solutions if Ki ^ 0, we also have a 
term that includes Ii, which is the conserved mass. We 
multiply ([38]) by /' and integrate again with respect to 
y to get the following nonlinear differential equation for 
the traveling waves: 

lf'j^^-{fTF = K^f + K,. (41) 

We see that K2 must also be zero for solutions / that 
are compact. Now, if we set Ki = and K2 = and 
integrate with respect to y, we obtain 



_ c 1 



1(1-1) 



(42) 



From ((M)) and ([i^ we can solve for Jjn,p and Ii in 
terms of l2- We obtain 



{l~2){m~ 1) 
2[p + m+ (m- \)V\ 



and 



_ l{l - l){p + im - 2)c 
'~ 2[p + m{l + l)-l] 



(43) 



(44) 
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Notice that when / ^ 2, c = 1 for consistency. This is 
related to the fact that the equation for u(x,t) becomes 
a hnear equation with propagation velocity c = 1. The 
energy of the solitary wave is given by 



H 



1 



1 



1(1^1) 



(45) 



and the momentum P ^ I2/2. From Eqs. (gS]) and ([iS]) 
we deduce that the energy, momentum, and velocity of 
the solitary wave are related by 



H = Pc/t 



(46) 



where r is given in (|22l 



A. Weak solutions 

We are interested in compacton solutions that are a 
combination of a compact function f{x) confined to a 
region (initially — xq < a; < and zero elsewhere). At 
the boundaries ±xo the function f{x) is assumed to be 
continuous but higher derivatives most likely are not. For 
there to be a weak solution we require that the jump in 



1(1-1) 



(fTF - Kif - K2 



(47) 



be zero when we cross from xq — e to xq + e. Since f{xo) 
is assumed to vanish, the requirement for a weak solution 
is 



Disc[{fTi^)nx)U - 0, 



(48) 



where Disc is the discontinuity across the boundary xq- 
This is always satisfied if there is no infinite jump in the 
derivative of the function. 

We are concerned mostly with the cases where the in- 
tegration constants are set equal to zero. We can then 
rewrite (l47l) as 



£ f 2-p _ _J_ 

2-' 1(1-1) 



Notice that for to = 2 we recover 



£ f2-p t 

2-' 1(1 



1 



(49) 



(50) 



which was studied previously. The solitary wave for (jSOp 
is obtained by joining the positive and negative solutions 
of the square root of and for m = 2n one again ob- 
tains a real solution by joining the positive and negative 
real parts of the solutions of the mth root of (|49|) ap- 
propriately shifted so that the maximum is at the origin 
2/ = 0. 

If we are looking for compactons, then the finitencss of 
the derivative when / — > requires that 



P<2, p<l. 



(51) 



The case to = 2 gives both compactons and real equa- 
tions. For TO = 2, when 



l^p + 2, 



(52) 



the width of the compacton is independent of the veloc- 
ity. In our previous discussion of scaling we found that 
when I = p + m, the width of the compacton is indeed 
independent of the velocity. 



B. Compacton solutions when m is an even integer 

Compacton solutions are constructed by patching a 
compact portion of a periodic solution that is zero at 
both ends to a solution that vanishes outside the com- 
pact region to give the weak solution described above. 
Let us look at the generalizations of the compacton equa- 
tion when we go from m ~ 2 to m = 4. Consider the case 
when p = 1 and let I = 3 and I ~ 4. For p = 1. Z = 3 the 
equation for the solitary wave is 



2-' 



(/')' 



For the positive branch of the solution, we get 

du 

la 



X — ct 



u - 



Performing the integral, we obtain 



ct = 2-3" 



-B 



'm— 1 m — 1> 



(53) 



(54) 



(55) 



where Bn(x,y) is the incomplete beta function while 
B(x, y) (see below) is the complete beta function. 
For TO = 2 this simplifies to 



x-ct = 2V6 



sm 



VfJW) 



which leads to the previous compacton result 



f = ic sin 



^(^-ct) 



For TO = 4 we get for the positive real fourth root 



x-ct^2''^'d^'^^c\B^(ll) 



(56) 



(57) 



(58) 



In Fig. [T] we plot B ^ 



'3 3^ 



- S(|,|) and its 



ror image as a function of f/(3c). Here, y ~ (x — 

ct)2-l/43-3/4^-l/2. 

Now consider the case p ^ 1 and I = 4, where the 
solitary-wave equation becomes 



For the positive branch of the solution, we get 

1 31-l/m 



X — ct — 



(59) 



(60) 
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We also find that 




f[y]/3c 



FIG. 1: f /{3c) versus y for m — 4, I — 3, and p = 1. 



or 



(61) 



By adding the positive-real- and negative-real-root solu- 
tions (for m an even integer) we get the complete com- 
pacton profile. The compacton vanishes elsewhere. 
For m ~ 2 this leads to 



x-ct=[3/(2c)]i/4B^(i,i) 



(62) 



When 111—2, the equation for the compacton can in fact 
be directly solved in terms of elliptic function. In partic- 
ular, consider the equation at m = 2 



2/ ^2-^ '■■^ ^ ■ 



(63) 



(64) 



Assuming a solution of the form 

f^Acn^(3y,e^l/2), 
we find that 

P = c^/\96)-^/\ A = V6^, (65) 

where cn(a;, k) is the Jacobi elliptic function with modu- 
lus k. 

C. Hyperelliptic compactons 

Consider next the generalization of the hyperelliptic 
compactons discussed in Ref. [4]. For this purpose we 
assume that we can parametrize the solutions by 

f = AZ-[l3{x - ct)], (66) 
and we demand that 

(Z')'" = l-Z^\ (67) 
This immediately leads to the relations 

m m(l — 2) 



m + p — 2 ' 



2(to-^P- 2)' 



(68) 



^ " ^ - 2^ ~ ^(^-1) 



This leads to 



[d(Z- l)/2]i/('-2) 



and 



a[^(^-l)]vJ -^(^-^)/^] 



(i_p_,„)/(„j(;_2)) 



(69) 



(70) 



(71) 



Note that this ansatz gives the correct scaling behavior 
of the amplitude parameter A and the width parameter 
P for the velocity c. 

The solution to the differential equation ((67)) has m 
branches corresponding to the various values of e^*'^"/™ 
when m is an integer and n = 1, 2, . . . , m. For even 
integer m, the positive root can be integrated to give 



/ dx ( 



1-x' 



- l/rn 



^ ^ ' m' 2r 2r' 



(72) 



where 2F1 denotes the hypergeometric function. For even 
m we get the full solution for the compacton by adding 
the positive-real-root and the negative-real-root solutions 
to get the complete compacton profile. The compacton 
vanishes elsewhere. 



V. CONSERVED QUANTITIES 



For solutions satisfying ((67|) it is possible to determine 
explicitly the conserved quantities in terms of the velocity 
c of the wave and the parameters and m of the 

differential equation. We have already found that the 
parameters A and /3 are given by (fTD]) and ((7T|) . 

We use the generic integral 



U j 2TV{^+l3 + iy ^^^^ 



For the mass we get 

M = J dxAZ^iPx) = ^ y dZidZ/dyr^Z" 

= 2- f dZZ°(l -Z^^)^!/™ 
P Jo 



A 

-4r(^)r(^) 



(74) 



Here, we have used the fact that the total area under the 
solitary wave is twice the area coming from the positive 



6 



root. (See Fig. [T]) For the momentum we get 



2/3rr (1 - + ' 

where a — m/(m + p — 2) and the energy is E — cP/r 



(75) 
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z[y] 



VI. SPECIAL CASES 

There are two types of special cases. The first occurs 
when T, I, and p are integers. From Eq. (|68p it follows 
that for m = 4 



I = 2 + t{p + 2)/2. 



(76) 



We consider here the set p = {1,2}. For p = 1, r = 2n 
and I = 3n + 2. For p = 2, we find instead that t, 
Z = 2n + 2, n = 1, 2, . . .. For m 6 



I = 2 + r(p + l)/3. 



(77) 



Thus, T = 3n, I = 2 + 7i(p + 1) yields simple integer 
solutions for integer n, p = I, 2. 

The other interesting case arises when the width of the 
solitary wave is independent of the velocity. This occurs 
when 



I =: p + m. 



(78) 



We study the cases for which {p,l} = {1,1 + m} and 
{p, 1} = {2,2 + to}. When Z = p + m, t is given by 



T = to/2 

and here we also have a = m/{l — 2). 



(79) 



Case m = 4 



For this case 



2/ -4 4 
T = — , a = 



(80) 



For T = 2, Z = p + 4 and the width is independent of 
velocity. The relevant function to invert is 



y 



dx 



(l-x4)l/4 



— Z 2F1 



,4' 4' 



(81) 



For a compacton centered about the origin y = 0, the 
two halves of the compacton are given by 



y±=±/i(Z)T/i(Z=l), 



(82) 



FIG. 2: Z[y] for m = 4, / = p + 4 r = 2. 



where 



h{x)^X2Fi{ 



and 



/i(^=l) 



r(l)r(|] 



i7^^/2 = 1.11072. 



(83) 



(84) 



The result for Z[y\ is shown in Fig. 2. 

For the case p = 1. Z = 5, we get A = (0/2)^/-^ and the 
solution goes as Z^/^ . For the case p = 2, / = 6, we get 
A = (c/2)^/^ and the solution is linear in Z. 

When r = 3; there is another special case with integer 
values: p ~ 2; I — %^ A — \. Then, the relevant function 
to invert is 



y = 
Since 



dxil 



„6^-l/4 



Z2F1 (i, i; 



(85) 



2-^ 1 V6' 4' 6' 



)-r(|)r(i)/r(ii), (86) 



we again use both the positive and negative solutions to 
make up the entire function Z(y). The case p = does 
not correspond to compacton solutions and p > 2 does 
not allow for a finite derivative when the amplitude of 
the solitary wave becomes zero. 



B. Case m = 6 

For integer r = 3n we discuss the r = 3 case. Here, 
I = p + 3, and there are two possibilities: Let p = 1, 
I = 4. Then, 



For p = 2, ; = 5, 



^cxc^/^ a = 6/5. 



A(xc^/^, a=l. 



(87) 



(88) 
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The form of the function Z[y] is now similar to the 
previous case: 

f2iZ) = r dx (1 - x')-'/' = Z,F,{l\-l-Z^). 

JO 

(89) 

We obtain the two halves of the function Ziy) by invert- 
ing 



y±^±h{Z)Tf2{Z^l), 



where 



/2(z = i) = r (I) r (I) = 1 = 1.0472. 



(90) 



(91) 



VII. ALTERNATIVE GENERATING 
FUNCTION AND STABILITY 

Solitary waves of the form f{y) — f{x — ct) can be 
derived by considering the following function: 



/'(?/)] = / dx{H[f,f]+P[f]c) 

dx^ifj']. (92) 



Notice that (p is the negative of the Lagrangian density. 
That is, the original equation for the solitary wave can 
be written as 



/5 n 



(93) 



The once-integrated equation ((38l) (with no integration 
constants) is obtained from 



6^ 



0, 



or equivalently from the Euler-Lagrange equation 
df d / dip 

df^d^ [W. 

We have explicitly 



dy 



(94) 



(95) 



(96) 



The first variation after an integration by parts can be 
written as 



5$ 



dy 



f 



(-1 



I - 1 



-mnrr-^f" 



Sf. 



(97) 



The second variation which is important for the linear 
stability analysis can be written as 



SH^ / dy6f L6f, 



(98) 



where L is the operator 

L = c - f-^ - Pip - l)fP-^{f'r - mpF-\fT-^f" 



{-mpr-\fT~' - ™(m - 2)r{fT'^f") 



dy 
(99) 



When m = 2, this reduces to the result given in Dey and 
Khare^. 



One can write $ in terms of /; and J„ 



m — 1 



~ Jrn . n 



1(1 



(100) 



Following Derricki^, we consider the scale transformation 
X — > Xx. Under this transformation 



1 



so that 
*[/(Ay)] 



1 



\m-l J 



II 



Xl{l - 1) 



2A 



(101) 



h. (102) 



If we assume that taking the derivative of (f> with respect 
to A and setting A = 1 gives a solution, we get 



d$(A) 



dX 



Ii 



(103) 



This is precisely the equation of motion integrated 
over space that we found earlier. [Derrick looked to see 
if the second derivative of (|103p became negative which 
would indicate that the solitary wave was unstable.] If 
we calculate the second derivative, we obtain 



d2$(A) 



c/2 



I ~ 2 



dX^ 21(1 - + m) + m + p] 



■n? -7m-2p + A)l + 6m + 2p~4]. (104) 

This does not factor to give a simple criterion for sta- 
bility. However, another choice leads to a simple stability 
criterion. Suppose we instead make the scaling 



f{y) ^ A''/(Ay). 



(105) 



This again leads to the equations of motion plus a bound- 
ary term because 



dA 



A=l 



dy 

dip 



df 



dip d dip 
df ^ diW 

{pf + x.f") 



" = 0. (106) 



Assuming that the boundary term vanishes at the edges 
of the compacton, we recover the equation of motion 



^XPf{Xy)] = 



\m— l+p(rn+p) t 

m — 1 

-^A'-^ + ^A^-^/.. (107) 
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The condition for a minimum is 



where 



dX 







TO — 1 + p(to + p) 
TO — 1 



-(2p-l)/2. (108) 



The particular case p = 1/2 is special in that the con- 
served momentum P is invariant under this transforma- 
tion; that is. 



P[Ai/2/(Ay)] = P[fiy)]. 



(109) 



If wc choose p ~ 1/2. when wc vary we are varying 
the Hamiltonian with the constraint that P is held fixed. 
This is exactly what happens in a trial variational calcu- 
lation where the parameter A, now thought of as a varia- 
tional parameter, is a constraint variable to be eliminated 
and determined in terms of the momentum P, which is 
the dynamical variable of the reduced Hamiltonian sys- 
tem. (This particular variation was first considered by 
Kuznetso\*ii and was then elaborated on by Karpman^ 
and Dey and Khare^.) 
For p= 1/2 we obtain 



(110) 



(/-2)(to-1) ^ 
;(?-l)(3m + p-2) 



which is precisely the relation obtained by using (j43|) and 
(im. For arbitrary p we get a linear combination of ([55]) 
and For arbitrary p the second derivative does not 

factor into a simple form that allows one to say when 
it changes sign. However, for p = 1/2 the answer does 
factor and the second derivative yields 



$"(A) 



Pc{l - 2){3m + p-l){'Sm+p~2) 



A=i 4;(-1-|-to)-|-to-|-p 

(111) 

We also learn from the conditions in (|?T|) that a weak 
solution that is compact can exist ii p < 2, p < I. This 
leads to the statement that solitary waves will be unsta- 
ble under this type of deformation when 



I > p + 3m. 



(112) 



More general scale transformations involving two param- 
eters, such as 

f{y) - /i'/'/(Ay), (113) 

have been discussed by Karpmann^ and Dey and Khare^. 

A. Linear Stability 

In this section we extend the analysis of Karpman^ to 
our generalized KdV equation. To study linear stability 
we assume that we can write 



if, a) 



dxf*g. 



(115) 



We change variables to y = x — vt and T = t, so that 
we parametrize an arbitrary addition to f{y) as v{y, T). 
The linearized equation of motion for v{y,t) is then 



±v{y,T)^-^[Lviy,T)], 



(116) 



where L is given by ((99)) . 

The equation l[37|) for the solitary wave can be written 

as 



Ldyfiy) = 0. 



(117) 



Thus, f'{y) is a zero eigenfunction of L corresponding 
to the translation invariance of the solution (that is, the 
Goldstone mode). 

Now, if we take the derivative of the first integral of the 
solitary wave equation ([551) with respect to the velocity 
c, we obtain 



L^--f 



Inverting, we get 



dl 
dc 



-L-\f. 



This result will be useful later. We also have 

(/, /')=o 



(118) 



(119) 



(120) 



from integrating by parts. If we now consider 

viy, T) = e-^-^^{y) + e^'^* V(2/), (121) 



then ip satisfies 



(122) 



When m is an even integer, L is a Hermitian operator. 
In that case there is a theorem that all the w are real if 
one of the two operators on the right side of p22p is pos- 
itive definitoi^. (The more general case will be studied 
elsewhere.) A sufficient condition for real eigenvalues is 



(123) 



where ip is orthogonal to / and /'. This condition is 
exactly the same as requiring that the second variation 
oi ^ = H + Pc be a minimum at /. If the solitary- wave 
solution is a minimum of <&, then the solution is linearly 
stable. 

Our objective is to find the extremal value of {ipi^Lip) 
and to find the criterion that guarantees that it is pos- 
itive. To find the extrema one solves the constrained 
variation condition 



u{x,t) = f{y) + v{x,t), |w|«l, (m,w)=0, (114) 



5[(V,i^)-A(V',V)-C^(V',/)] =0, (124) 
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so that 

(L - A)^ = Cf. 
Using {ip, /) ~ 0, wc find that 



(125) 



(126) 



One solves this equation by expanding ip and / in a 
series of eigcnfunctions of the operator L. Letting L0„ = 
A„0„ and assuming the ordering > Am if n > m, we 
find that Ai = and 0i = /'. Letting 



/ ^ bn<l>n, 



(127) 



we find that 



From (-0, /) = we obtain 



A„ - A 



(128) 



(129) 



We are interested in the lowest solution A,„i„ that solves 
this equation because this gives the minimum of {tp,Lip) 
that we seek. Assuming along with Karpman^ as well as 
Dey and Kharo^ that L has only one negative eigenstate, 
one then deduces that Lip) > is satisfied if r(0) < 0. 
However, when A we get 



L^A=0 = Cf. 

Thus, from pTH]) we find that 

df 



dc 



and 



This criterion, namely that 



dP 



dP 
dc 



(130) 



(131) 



(132) 



(133) 



for stability, gives exactly the same sufficient result for 
stability as all the other criteria used. 

Since for all of our solutions P oc c(p+3"-0/[('-2)™] ^ 
it immediately follows that these solutions are linearly 
stable provided that 



2 <l <p + 3m. 



(134) 



B. Lyapunov Stability 

Lyapunov stability uses sharp estimates and has been 
used by Weinsteini^ and Karpman et al.— . Here we 
want to show that the compacton solution is a mini- 
mum of the Hamiltonian for fixed momentum P. We 
show this using Holder's Inequalitj^ and follow the ar- 
guments of Weinsteirti^, KuznetzoAiii, Karpmann^, and 
Dey and Khare^. We do this by first showing that the 
Hamiltonian for fixed momentum P is bounded below 
and then that the compacton solution satisfies the con- 
dition that it is a particular lower bound. We can write 
the Hamiltonian H in terms of /; and Jm,p as 



H[f] 



where 



m — 1 



h = 2P. 



1(1^1) 



Now consider that 



// < (max / 



l-2\ 



(135) 



(136) 



(137) 



We want to bound /; by a function of Jm.p and the con- 
served momentum P. To do this it is convenient to write 



jl-2 ^ jjaj(i-2)/a 

Then, writing 

r = j dydr/dy = a J dy f 
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(138) 



a-l-k ^/ rk 



dy 



r, (139) 



we use the Holder incqualityii to show that 



r < a(^| dy\f 



{a—l — k)rn^jf'yrL 



l/m 



X ( / dy\f^{y)\ 



(140) 



We can dispense with the absolute values when m is an 
even integer. Thus by choosing 



kj 



p = m{a — 1 — fc), 



(141) 



we can relate the second term of (|135p to P and Jm,p- 

Moreover, to relate the bound to the energy of the 
solitary wave (rather than having a general lower bound 
that depends on the choice of a) , one must further choose 



ma = p + 3m — 2, 
from which we find that 



771—1 



(142) 



(143) 
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Doing this we obtain 



The stationarity condition is 



H > niin 



m — 1 



where 



and 



(/ - 2)(m- 1) + P + 3to- 2 
p + 3m — 2 ' 



1 + 3m - 2 



l{l - 1) V m 



m(;-2)/(p+3in-2) 



(144) 
(145) 

(146) 



Minimizing with respect to Jm,p for fixed momentum P, 
we obtain 



I — p — 3m 
(m- 1)(/ - 2) 



(147) 



For the sohtary wave which obeys the generahzed KdV 
equation we get 



J, 



{l-2){m- 1) 
p + m + (m — 1)1 



PsolC 



(148) 



Thus, the sohtary wave is a minimum of the Hamihonian 
and satisfies 



Esol — Ejnin — PsolC / 



(149) 



where r is as given by (1^^ . Thus, as long as 2 < Z < 
p + 3m, one has stable solitary waves. 



VIII. VARIATIONAL STABILITY OF 
SOLUTIONS 

Suppose that we have found an exact solution of the 
form AZ{f3[x — q{t)]}. Then one can find sufficient con- 
ditions for instability of this type of solution by seeing 
if the solution is a minimum rather than a maximum of 
the Hamiltonian as a function of /? with the conserved 
momentum P held fixed. 

We know that the exact solutions are stationary for 
fixed P under variations in /3. That is. 



dH 
'dp 



= 0. 



(150) 



We can write the Hamiltonian in the generic form (P 
fixed) 



(151) 



where the constants depend on P and the parameters 
that define the Lagrangian and a and b also depend on 
the parameters in the Lagrangian. 



f)M 1 

^ = - - [-C^aP^ + C2b/3''] , (152) 

from which we infer that 

CiaP" = C26/3^ (153) 
The edge of stability of these solutions is given by 



(9/32 



= = -Cia(a- l)/3° + C26(6- l)/3^ (154) 



At the minimum this leads to the condition that 

a^{l-2)/2, 6=(p + 3m-2). (155) 
Thus, the critical case is 

I ^ p + 3m, (156) 



which agrees with (|134p . We expect that the solutions 
we have found are stable as long as 



I < p + 3m. 



A. Approximate variational solutions 



(157) 



To study stability it is useful to have approximate so- 
lutions that are close to the exact solutions to see if they 
relax to the exact solutions or become unstable. For this 
purpose it is useful to study the post-Gaussian trial func- 
tions 

fvix - ct) EE g{x - ct) = Aexp [~\l3{x - ci)|^"] , (158) 

where /?, and n are continuous variational parameters 
chosen to minimize the action. These trial wave func- 
tions have earlier been successfully usedi^ to approximate 
various solitary waves in both KdV systems and NLSE 
applications. 

The advantage of these trial functions is that the action 
as well as all of the conserved quantities can be explicitly 
evaluated using the formula 



l/2n' 



(159) 



In terms of our previous notation we have 

Z(z) ^cxp(-|z|2") . (160) 
Thus, using the trial function, we get 



Ci(0 = ldzZ\z)^2j dze-'l^l' 

= r^r(i + ;l 

2n 



(161) 
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We also have 

C2(p, m) 



dz [Z'{z)]"'ZP 



2 / dz(-2)™e- 
-(-2n)'"(m+p) 



mx ^ —p ^m,{2n—l) 



'T (m 



2n I 



.(162) 



These expressions depend on the variational parameter 
n, which determines the shape of the solution; n = 1/2 
gives a peakon shape and n = 1 is the usual Gaussian. 
From these we can determine the quantities 



dip, to) = K{p, TO, n)T (1 + 
where 



r TO 



m-l \ 
2n ) ' 

(163) 



+ m + 2np+p 



Kij), TO, n) - 
and 



(to+p) 2„ 



TO — 1 



^"^'-^--r(i + ^)- 



/-I 

For the trial function we get 



and 



V + m-l 

(3 = Pi-p-s^ 



Ciip, m){p + 3to — 2) 
^3(0(^-2) 



(164) 
(165) 

(166) 
(167) 



To determine the best trial function in this class we must 
also minimize the Hamiltonian with respect to the pa- 
rameter n. As in our discussion of conservation laws, we 
get 



H = f{l,p,m,n)p-\ 



(168) 



where r is given in (22) and f{l,p,m,n) is given in ([S^) - 
For solutions that are compact and cover half of the pe- 
riod of a positive periodic function, an alternative choice 
for a variational trial function is 



U2{x) = A[cos{l3x)P , 



(169) 



where (3 and 7 are the variational parameters. For inte- 
ger p, TO, and / it is again possible to obtain an explicit 
expression for H [(],"/]. One can perform the minimiza- 
tion with respect to [3 explicitly. Determining the global 
minimum in the parameter 7 must be done numerically. 



B. Case p = 1, I = 3, m 

First, consider to = 2, where the exact solitary wave 
solution is 



f{y) = 3csin^ 



y 

276 ^2 



(170) 




FIG. 3: f{y) versus y for m — 2, I — 3, p — 1. This graph 
compares the exact versus variational solitary waves for P = 
1. 



For simplicity we normalize our functions by choosing 
P = 1. Then, this solution is 



/(Z) = 25/^3-3/^-1/2 ^os^ 



2^/6 



(171) 



This belongs to the class of variational solutions of the 
second type, and we would have obtained this exact an- 
swer from our variational minimization procedure. For 
the post-Gaussian trial functions, the lowest-energy vari- 
ational solution having P = 1 is found to be 



u{x) = 0.583578 exp (-0.0314705 x 



2.308\ 



(172) 



Fig. 3 compares the exact and variational function for 
P = 1. Note that apart from the region where the com- 
pacton goes to zero the agreement is excellent. 

For TO = 4 the value of n that minimizes the first trial 
function is n = 0.920655. Again, normalizing to P = 1 
(which yields c « 1/3), we find that the best hmction in 
this class is: 



u{x) ^ 0.995936 exp (-0.396108 a;^ ®''"^) 



(173) 



If we use the second type of trial function we find that 
the values of P and 7 that give a global minimum in the 
reduced space are 



(3 = 0.342787, 7 = 5.67846, 



(174) 



which leads to ^ = 0.97067 for P = 1. Thus, the best 
trial function in the second class is given by 



U2ix) = 0.97067 [cos(0. 3427872;)] 



5.67846 



(175) 



In Fig. 4 we compare the two variational approxima- 
tions and note that apart from the fact that the u{x) is 
not compact, the agreement is quite good. Both solutions 
are global minima of the respective reduced Hamiltoni- 
ans, which depend on two parameters. 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 

FIG. 4: u{x) and U2{x) versus x for m = 4, I = 3, p = 1. 




FIG. 5: f{y)/3c versus y for m — 4, I = 3, p — 1. Exact 
answer versus two variational solitary wave solutions. 



To compare our implicit exact result with our vari- 
ational approximations, we change variables to y = 
x/2^/^3'^/^c^/^ and redisplay the approximate solitary 
waves shown in Fig. 5 along with the exact result of 
Fig. 1. 
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